Noncolliding Brownian motions and 
Harish-Chandra formula 



Makoto Katori 
Department of Physics, 
Faculty of Science and Engineering, 
Chuo University, 
Kasuga, Bunkyo-ku, 
Tokyo 112-8551, Japan 
e-mail: katori@phys.chuo-u.ac.jp 

HlDEKI TANEMURA 

Department of Mathematics and Informatics, 
Faculty of Science, 
Chiba University, 
1-33 Yayoi-cho, Inage-ku, 
Chiba 263-8522, Japan 
e-mail: t anemur a@ mat h . s . chiba-u . ac . j p 

Abstract. We consider a system of noncolliding Brownian motions introduced in our pre- 
vious paper, in which the noncolliding condition is imposed in a finite time interval (0, T\. 
This is a temporally inhomogeneous diffusion process whose transition probability density 
depends on a value of T, and in the limit T —>■ oo it converges to a temporally homoge- 
neous diffusion process called Dyson's model of Brownian motions. It is known that the 
distribution of particle positions in Dyson's model coincides with that of eigenvalues of a 
Hermitian matrix-valued process, whose entries are independent Brownian motions. In the 
present paper we construct such a Hermitian matrix-valued process, whose entries are sums 
of Brownian motions and Brownian bridges given independently of each other, that its eigen- 
values are identically distributed with the particle positions of our temporally inhomogeneous 
system of noncolliding Brownian motions. As a corollary of this identification we derive the 
Harish-Chandra formula for an integral over the unitary group. 
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1 Introduction 

Dyson introduced a Hermitian matrix- valued process whose r/'-entry equals to By (t) / \/2 + 
y/—lBij(t)/\/2, iil<i<j<N, and equals to Bu(t), if i = j, where B i: j(t), B^ (t), 
1 < ^ < j < iV, are independent Brownian motions [3]. He found that its eigenvalues 
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perform the Brownian motions with the drift terms acting as the repulsive two-body forces 
proportional to the inverse of distances between them, which is now called Dyson's model 
of Brownian motions. A number of processes of eigenvalues have been studied for random 
matrices by Bru 00], Grabiner [Ij, Konig and O'Connell [T%] . and others, but all of them are 
temporally homogeneous diffusion processes. In the present paper we introduce a Hermitian 
matrix- valued process, whose eigenvalues give a temporally inhomogeneous diffusion process. 

Let Y(t) = (Yi(t), Y 2 (t), . . . ,Y N (t)) be the system of N independent Brownian motions 
in R conditioned never to collide to each other. It is constructed by the /i-transform, in the 
sense of Doob [I], of the absorbing Brownian motion in a Weyl chamber of type Ajv-i, 

R^ = {x G R N : x x < x 2 < ... < x N } (1.1) 

with its harmonic function 

l<i<j<N 

x G R< . We can prove that it is identically distributed with Dyson's model of Brownian 
motion. In our previous papers we introduce another system of noncolliding Brow- 

nian motions X(t) = (Jfi(t), -X^t), . . . , Xj^(t)), in which Brownian motions do not collide 
with each other in a finite time interval (0, T]. This is a temporally inhomogeneous dif- 
fusion process, whose transition probability density depends on the value of T. It is easy 
to see that it converges to the process Y(£) in the limit T — > oo. Moreover, it was shown 
that P(X(-) G dw) is absolutely continuous with respect to P(Y(-) G dw) and that, in 
the case X(0) = Y(0) = 0, the Radon-Nikodym density is given by a constant multiple of 
l/hj^(w(T)). Since this fact can be regard as an iV-dimensional generalization of the relation 
proved by Imhof jU] between a Brownian meander, which is temporally inhomogeneous, and 
a three-dimensional Bessel process, we called it generalized Imhof s relation jToj. 

The problem we consider in the present paper is to determine a matrix-valued process 
that realizes X(t) as the process of its eigenvalues. We found a hint in Yor [21] to solve 
this problem : equivalence in distribution between the square of the Brownian meander and 
the sum of the squares of a two-dimensional Bessel process and of an independent Brownian 
bridge. We prepare independent Brownian bridges /3ij(t), l<i<j<Nof duration T, which 
are independent of the Brownian motions Bij(t), 1 < i < j < N, and set a Hermitian matrix- 
valued process S T (t), t G [0,T], such that its ij-entry equals to Bij(t)/\/2 + y/^-iPij(t)/^/2, 
if 1 < i < j < N, and it equals to B„(t), if i = j. Then we can prove that the eigenvalues 
of the matrix S T (t) realize X(t),t G [0, T] (Theorem 2.2). This result demonstrates the fact 
that a temporally inhomogeneous diffusion process X(t) in the iV dimensional space can be 
represented as a projection of a combination of iV(iV + l)/2 independent Brownian motions 
and iV(iV — l)/2 independent Brownian bridges. 

It is known that Brownian motions Bij(t), 1 < i, j < N are decomposed orthogonally into 
the Brownian bridges B^t) — (t/T)Bij(T) and the processes (t/T)Bij(T) (see, for instance, 
j2Hl[24j). Then the process S T (t) can be decomposed into two independent matrix- valued 
processes Q^\t) and Q^ 2 \t) such that, for each t, the former realizes the distribution of 
Gaussian unitary ensemble (GUE) of complex Hermitian random matrices and the latter 
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does of the Gaussian orthogonal ensemble (GOE) of real symmetric random matrices, re- 
spectively. This implies that the process S T (t) is identified with a two- matrix model studied 
by Pandey and Mehta j20J |22|; which is a one-parameter interpolation of GUE and GOE, 
if the parameter of the model is appropriately related with time t. In [IA we showed this 
equivalence by using the Harish-Chandra formula for an integral over the unitary group [8 u . 
The proof of Theorem 2.2 makes effective use of our generalized version of Imhof's rela- 
tion and this equivalence is established. The Harish-Chandra formula is then derived as a 
corollary of our theorem (Corollary 2.3). 

As clarified by this paper, the Harish-Chandra integral formula implies the equivalence 
between temporally inhomogeneous systems of Brownian particles and multi-matrix models. 
This equivalence is very useful to calculate time-correlation functions of the particle systems. 
By using the method of orthogonal polynomials developed in the random matrix theory , 
determinantal expressions are derived for the correlations and by studying their asymptotic 
behaviors, infinite particle limits can be determined as reported in |2Tj IT3] . 

Extensions of the present results for the systematic study of relations between noncol- 
liding Brownian motions with geometrical restrictions (e.g. with an absorbing wall at the 
origin [T3| I17|) and other random matrix ensembles than GUE and GOE (see |19| 123} IT], for 
instance), will be reported elsewhere [TB] . 



2 Preliminaries and Statement of Results 
2.1 Noncolliding Brownian motions 

We consider the Weyl chamber of type A^_i as [51 [7j. By virtue of the Karlin-McGregor 
formula fJJE], the transition density function fw(t, y|x) of the absorbing Brownian motion 
in R< and the probability A/jv(t, x) that the Brownian motion started at x e R< does not 
hit the boundary of R< up to time t > are given by 



f N (t,y\x) = det G t ( Xj ,yi) , x,y e R£, (2.1) 

l<ij<iV 



and 



J\f N (t,x.) = / dyf N {t,y\x) 



respectively, where G t (x,y) = (2nt)~ 1 ^ 2 e~^ y ~ x ^ l 2t . The function h^(x) given by (jl.2j) is a 
strictly positive harmonic function for absorbing Brownian motion in the Weyl chamber. We 
denote by Y(t) = (Yi(t),Y 2 (t), . . . ,Y^(t)),t G [0, oo) the corresponding Doob /i-transform 
|lj , that is the temporally homogeneous diffusion process with transition probability density 
p N (s,x,t,y); 



£-N 2 /2 



y 



p N (0,0,t,y) = — — exp -^f } h N (y) 2 , (2.2) 



d(N) ^ { 2t 
1 

h N (x 



Pn(s, x, t, y) = f N (t - s, y\x)h N (y), (2.3) 
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for < s < t < 00, x, y G R<, where Ci(N) = {2<k) n I 2 ]jf =1 T(j). The process Y(t) 
represents the system of iV Brownian motions conditioned never to collide. The diffusion 
process Y(i) solves the equation of Dyson's Brownian motion model [3] : 

dYi{t) = dBi(t) + vm-Y-ti) ^ tG [°'°°)' * = 1, 2, (2.4) 

where Bi(t), i = 1, 2, . . . , N, are independent one dimensional Brownian motions. 
For a given T > 0, we define 

A/iv(T - s,x) 

for < s < i < T, x, y G R<. The function g%(s, x, t, y) can be regarded as the 
transition probability density from the state x G R< at time s to the state y G R< 
at time t, and associated with the temporally inhomogeneous diffusion process, X(£) = 
(X%(t), X 2 (t), . . . , X N (t)), t G [0,T], which represents the system of N Brownian motions 
conditioned not to collide with each other in a finite time interval [0, T\. It was shown in 
[To] that as |x| — > 0, g^(0,x, t,y) converges to 

T N(N-l)/4 t -N 2 /2 f \ \2> 

gl(0, 0, t, y) = —- exp |-J|i- J h N (y)Af N (T - t, y), (2.6) 

where C 2 (iV) = 2 W / 2 [Xjii r(j/2). Then the diffusion process X(£) solves the following 
equation: 

dXi(t) =dBi(t) + bJ(t,X(t))dt, te[0,T], i = 1, 2, . . . , JV, (2.7) 

where 

6f(*, x) = — In A^(T - t, x), i = 1, 2, . . . , N. 

From the transition probability densities (|2.2)1 . (J2.H)) and (|2.fij) . (|2.5jl of the processes, we 
have the following relation between the processes X(t) and Y(i) in the case X(0) = Y(0) = 

P X • e d«; = ^ L (t 7 rrvT p Y e A* . (2.8) 
C 2 (N)h N (w(T)) 

This is the generalized form of the relation obtained by Imhof [5] for the Brownian meander 
and the three-dimensional Bessel process. Then, we call it generalized Imhof s relation. 

2.2 Hermitian matrix- valued processes 

We denote by TC(N) the set of N x N complex Hermitian matrices and by S(N) the set of 
N x N real symmetric matrices. We consider complex- valued processes Xij(t), 1 <i,j < M 
with Xij{t) = Xji{t)\ and Hermitian matrix-valued processes S(i) = (^y(t))i<i,j<jv- 
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Here we give two examples of Hermitian matrix- valued process. Let B^(t), Bjj(t), 1 < 
i < j < N be independent one dimensional Brownian motions. Put 

4yflg(f), ifi<j, f 4=4(*)» if *"<i> 

xf j (t)={ ^2 and ^ (t)= J V2 

BfM, ifi = i, 1 0, ifi = j, 

with xfj(t) = xfi(t) and x\j(t) = —x l ^{t) for z > j. 



(i) GUE type matrix-valued process. Let H GUE (t) = (xf^t) + \/-T^L(0)i<tj'<JV- For 
fixed t G [0, oo), S GUE (t) is in the Gaussian unitary ensemble (GUE), that is, its probability 
density function with respect to the volume element U(dH) of 7i(N) is given by 

f -N 2 /2 f i 1 

/i GUE (#, t) = —— exp <^ --TrH 2 }, He H(N), 



where C 3 (N) = 2 n ^ 2 tt n2 ^ 2 . Let U(iV) be the space of all AT x AT unitary matrices. For 
any U G U(AT), the probability fi GlJE (H,t)U(dH) is invariant under the automorphism 
H — > U^HU . It is known that the distribution of eigenvalues x G R< of the matrix ensembles 
is given as 

r(x,t)~exp 4 Mx) 2 , 



Cx(iV) 2t 
dU, and so p^(0,0,t,x) = # GUE (x,t) from Q . 

(ii) GOE type matrix- valued process. Let S GOE (t) = (^(t))i<ij<jv- For fixed £ G 
[0, oo), S GOE (t) is in the Gaussian orthogonal ensemble (GOE), that is, its probability density 
function with respect to the volume element V(dA) of S(N) is given by 

V GOE (At)= — W) exp|--TM 2 |, AeS(N), 

where C±(N) = 2 N / 2 7r N ( N+1 ^ 4 . Let be the space of all iV x A" real orthogonal matrices. 
For any V G 0, the probability fi GOE (H, t)V(dA) is invariant under the automorphism 
A — ► V T AV. It is known that the probability density of eigenvalues x G R< of the matrix 
ensemble is given as 

PSj, and so ^(0,0, t,x) = c/ GOE (x,£) from p^jl . 

We denote by U(t) = (uij(t))i<ij<^ the family of unitary matrices which diagonalize 
S(f): 

U(tf~(t)U(t) = A(t) = diag{A*(t)}, 

where {Aj(t)} are eigenvalues of E(t) such that Xi(t) < A 2 (t) < • • - < Ajv(t). By a slight 
modification of Theorem 1 in Bru [2] we have the following. 
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Proposition 2.1 Let Xij(t), 1 < i,j < N be continuous semimartingales. The process 
\{t) = (Ai(t), X 2 (t), . . . , X^it)) satisfies 

d\ i (t) = dM i {t) + dJ i {t), i = l,2,...,N, (2.9) 

where Mj(t) is the martingale with quadratic variation (Mi) t = J^Tu(s)ds, and Ji(t) is the 
process with finite variation given by 

N 1 

dJi ® = JL K (t)- Xj (t) 1{Xi ^ Xi)Tij{t)dt 

+ the finite variation part of (U(tydS(t)U(t))a 

with 

T l3 (t)dt = (U\t)d~(t)U(t)) l3 (U\t)d^(t)U(t)) 3i . (2.10) 

For the process 2 GUE (t), rfSg UE (t)dS^ UE (t) = 5 ie 5 jk dt and r y (t) = 1. The equation (Q 
is given as 

dXiit) = dB.it) + V 1 dt, 1 < i < AT. 
rrf. <M J J — A,(tJ 

Hence, the process X(t) is the homogeneous diffusion that coincides with the system of 
noncolliding Brownian motions Y(i) with Y(0) = 0. 

For the process H GOE (t), c/Sg OE (t)rfH^ OE (t) = \(6u8jk + 8ik8ji)dt and T i:j (t) = §(l + <%). 

The equation ()2.9|) is given as 

dMt) = dftW + \ £ A[(t) ! A (t) ^. i < i < Jv. 

2.3 Results 

Let /3ij(t), l<i<j<Nbe independent one dimensional Brownian bridges of duration T, 
which are the solutions of the following equation: 

M) = B\ 3 {t) - jf ^±ds, 0<t<T. 



For t G [0, T], we put 

1 



^M), if * < J, 
0, if i = j, 



with = — for i > j. We introduce the W(iV)-valued process H (t) = (x^-(t) + 

\/— l£y (t) ) x<i,j<N ■ Then, the main result of this paper is the following theorem. 
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Theorem 2.2 Let \i(t), i = 1,2, ...,N be the eigenvalues ofE T (t) with \i(t) < A 2 (t) < 
••• < Ajv(t). The process A(i) = (Ai(i), A 2 (£), . . . , Ajv(t)) i/ie temporally inhomogeneous 
diffusion that coincides with the noncolliding Brownian motions X(t) X(0) = 0. 

As a corollary of the above result, we have the following formula, which is called the 
Harish-Chandra integral formula jH] (see also JOJEI]). Let dU be the Haar measure of the 
space U(iV) normalized as fjj, N \ dU = 1. 

Corollary 2.3 Let x = (xx, x 2 , . . . ,x N ), y = (y x , y 2 , . . . , Vn) G R< • JVien 

/ dt/ exp (--Ltt(A x - tfA y U) 2 ) = f 1 ^" , det \G^ix uyj ) 

where A x = diag{xi, . . . ,xn} and A y = diag{|/i, . . . , 2/jv}. 
Remark Applying Proposition 12. II we derive the following equation: 



iX.it> = dB.lt) + -£ it L2__ *, (2.11) 



i = 1, 2, . . . , N, where U(t) is one of the families of unitary matrices which diagonalize S T (t). 
From the equations (J2.7)) and (j2.11j) we have 



GOE 



{dA){U{tfAU{t))i 



S(N) 

{1-M N {T -t,X(t)) ™ l ] 

- W + < r " " ( \( r - t ,A W ) " S^Fl) ■ (2 ' 12) 

The function A/jv(t, x) is expressed by a Pfaffian of the matrix whose r/'-entry is ^{{xj — 
Xi)/2y/i) with ty(u) = Jq e~ v dv. (See Lemma 2.1 in |15j.) Then the right hand side of 
(|2.12jl can be written explicitly. 



3 Proofs 

3.1 Proof of Theorem EH 

For y e R and 1 < i,j < N, let (3fj(t) = /3?.(t : y), t E [0, T], (j = R, I, be diffusion processes 
which satisfy the following stochastic differential equations: 

h h f* /$,•(« : y) - y 
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These processes are Brownian bridges of duration T starting form and ending at y. For 
H = (y§ + v / ^T^)i< 4J <iv e H{N) we put 



±=/3$(t:y/2y%), if z < j, 



&-yl), if i = ^, 

■ -^/%(*:^), ifi<J, 

0, if i = j, 

with £g(t : yfj) = Q(t : yfl) and : yL) = : y}<) for i > j. We introduce the 

^(iV)-valued process S T (t : H) = ($(t : yg) + V-lgfi : 4))i< ij <iv, t e [0,T]. From the 
equation (j3.1j) we have the equality 

E T (t : if) = ~ GUE (t) - f ET{S - H) - H ds, t G [0,T]. (3.2) 

Jo T -s 

Let Hu be a random matrix with distribution yU GUE (-,T), and Ao be a random matrix with 
distribution /i GOE (-,T). Note that /^(t : F), t G [0, T] is a Brownian motion when F is a 
Gaussian random variable with variance T, which is independent of B^(t),t G [0, T\. Then 
when ii[7 and A Q are independent of E GVE (t),t G [0, T], 

H T (t:fl £r )=S auB (0, «e[0,n (3.3) 
S r (t : A ) = ~ T (t), tG[0,T], (3.4) 

in the sense of distribution. Since the distribution of the process 5 GUE (t) is invariant under 
any unitary transformation, we obtain the following lemma from ()3.2|) . 

Lemma 3.1 For any U G U(iV) we have 

U^E T (t : H)U = E T (t :WHU), £g[0,T], 

in distribution. 

From the above lemma it is obvious that if H^ l > and are N x N Hermitian matrices 
having the same eigenvalues, the processes of eigenvalues of E T (t : H^'),t G [0, T] and 
E(t : H^),t G [0,T] are identical in distribution. For an N x N Hermitian matrix H 
with eigenvalues {ai}i<j<7v, we denote the probability distribution of the process of the 
eigenvalues of E T (t : H) by Qo a (-),t G [0, T]. We also denote by Q GUE (-) the distribution of 
the process of eigenvalues of E GVE (t),t G [0, T], and by Q T {-) that of E T (t),t G [0, T]. From 
the equalities (j3.3j) and ()3.4j) we have 



Q GUE (")= / Q£ a (-)<7 GUE (a,T)da, 
Q T (-)= / QU)g GOE (*,T) da . 
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Since Q GUE (-) is the distribution of the temporally homogeneous diffusion process Y(i) which 
describes noncolliding Brownian motions, by our generalized Imhof's relation ()2.8|) we can 
conclude that Q T (-) is the distribution of the temporally inhomogeneous diffusion process 
X(£) which describes our noncolliding Brownian motions. | 



3.2 Proof of Corollary 1231 

By (12.6)) we have 

$(0,0, t,y) 



C 2 (N) 
dz det 

R iV l<i,j<N 



1 rpN(N-l)/4:£-N 2 /2 



exp 



2?r(T - i) 



exp 



2* 

2(T-t) 



My) 



1 T N(N-i) /H -N>/2 {27r{T _ t ))-Wh N (y) 



dz det 



exp < — 



2t 2(T - 1) 



1 ^-i)/^/> (MT _ t))-^ My) 



dz exp 



, det 

2T | i<ij<iv 



exp 



T 



2t(T - 1) 



y.i 



T 



Setting (t/T)zi = a iy i = 1, 2, . . . , N, t(T — t)/T = a 2 and T/t 2 = a, we have 



$(0,0, t,y 



(2vr 



.-N/2 



N a N(N + l)/4 h 



X 



C 2 (N) 

/ (iaexp { |a| r det 

J k n I 2 J l<i,j< 



exp 



2a 



(3.5) 



We write the transition probability density of the process S T (t) by q%(s, Hi, t, H 2 ), 
< s < t < T, for H X ,H 2 E H(N). Then by Theorem E2 and the fact that U{dH) = 
Cu(N)h N (y) 2 dUdy, with C V (N) = C 3 (N)/d(N), we have 



g N (0,0,t,y) = Cu(N)h N (y) 



V(N) 



dU q%(0,O,t,lPA y U), 



(3.6) 



where O is the zero matrix. We introduce the H(N)-valued process Q^(t) = (of, ? (t))i<i,j<N 



and the 5(iV)-valued process Q^{t) = {dfj 1 {t))i<i,j<N which are defined by 



V2 



^1 



— \Bf 3 {t)--Bf ] {T)\ + ^_(3 lJ {t), ifi<j 



V2 



t 



B$(t) ~ j;B*(T), 



if i = j, 



and 



V2T 



5g(T), if i<j, 



respectively. Then E T (t) = 0«(t) + B^(t). Note that Bg(t) - (t/T)B%(T) are Brownian 
bridges of duration T which are independent of (t/T)B^(T). Hence Q^\t) is in the GUE and 

&W{t) is in the GOE independent of 9«(t). Since ^^(t) 2 ] = a 2 and ^[^(i) 2 ] = 1/a, 
the transition probability density q%(0, O, t, H) can be written by 



<&(0,O,t,H) 



V{dA) li 



GOE 



S(N) 



A,- )n 

a 



GUE, 



H — A,a 



C (iV)a- jv2 a jV ( jV + 1 )/ 4 
C 3 {N)C 4 (N) 



da h N (a) exp <! -^|a| 2 - — 
R iv I 2 Za z 



Tr(#-A a )H, (3.7) 



where we used the fact V(dA) = Co{N)hN{&)dVd& with the Haar measure dV of the space 
normalized as J dV = 1, and C (iV) = C^{N)/C 2 {N). Combining (jR~Hjl and (pT7jl 

we have 



d(N)<7 



N 2 -N 



(2n) N /*h N (y 

da/ijv(a) exp 



' da exp { |a| 2 ]" det 

Riv L 2 J l<i,j<N 



exp 



2a 2 



} /r jv rft/6XP {~2^ Tr(f/tAyf/ " Aa)2 } ' (3 ' 8) 



For each a > 0, ()3.8|) holds for any a > and we have 



A' 2 



C ^ det 
h N (y)h N (a) i<i,j<N 



V2 



: exp 



7T(J Z 



2<X 2 



- a i) 2 



U(JV) 



df/exp |--LTr(C/tA y [/ - A a ) 2 | 



This completes the proof. 
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